%% dx = (0.125,0.125)
solution = [
1 1.13315 1.28403 1.45499 1.64872 1.86825 2.117 2.39888 2.71828 
1.13278 1.28379 1.45489 1.64861 1.86753 2.11566 2.39754 2.71708 3.07922 
1.2807 1.45164 1.64542 0 0 0 2.71012 3.07167 3.48129 
1.44235 1.63499 0 0 0 0 0 3.45936 3.92071 
1.61515 1.83056 0 0 0 0 0 3.874 4.39042 
1.79517 2.03468 0 0 0 0 0 4.30578 4.87977 
1.97712 2.24081 2.53974 0 0 0 4.18562 4.7427 5.37436 
2.15447 2.44158 2.76693 3.13554 3.55268 4.02573 4.5612 5.16831 5.85645 
2.31978 2.62865 2.97865 3.37526 3.82467 4.33391 4.91097 5.56486 6.30581 
];
exact = [
1 1.13315 1.28403 1.45499 1.64872 1.86825 2.117 2.39888 2.71828 
1.13278 1.28361 1.45452 1.64819 1.86764 2.11631 2.3981 2.7174 3.07922 
1.2807 1.45122 1.64445 0 0 0 2.71123 3.07223 3.48129 
1.44235 1.63439 0 0 0 0 0 3.46001 3.92071 
1.61515 1.8302 0 0 0 0 0 3.87453 4.39042 
1.79517 2.03419 0 0 0 0 0 4.30638 4.87977 
1.97712 2.24036 2.53867 0 0 0 4.18555 4.74285 5.37436 
2.15447 2.44133 2.76639 3.13473 3.55212 4.02507 4.56101 5.1683 5.85645 
2.31978 2.62865 2.97865 3.37526 3.82467 4.33391 4.91097 5.56486 6.30581 
];
err = [
0 0 -4.44089e-16 -3.33067e-15 1.22125e-14 -1.84297e-14 -3.41949e-14 1.42109e-14 0 
-1.9984e-15 0.000181526 0.000375658 0.000421153 -0.000111815 -0.000654706 -0.000550191 -0.000313368 0 
4.44089e-16 0.000416533 0.000974841 0 0 0 -0.00110922 -0.000563375 0 
2.66454e-15 0.000593285 0 0 0 0 0 -0.000654101 -2.66454e-15 
-2.44249e-15 0.000354729 0 0 0 0 0 -0.000533956 8.88178e-16 
1.13243e-14 0.00049498 0 0 0 0 0 -0.000597157 4.44089e-15 
2.22045e-16 0.000443957 0.00107727 0 0 0 6.65158e-05 -0.000152307 0 
-2.70894e-14 0.000247099 0.000535565 0.00080784 0.000564554 0.000654141 0.000187694 1.35431e-05 0 
0 0 0 7.10543e-15 7.37188e-14 -6.92779e-14 2.66454e-14 0 0 
];
